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{^jQ Abstract 

r^ We study two-component ferniionic few-body systems in a one dimensional harmonic trap. Dif- 

^ ferent component fermions interact via a short-range interaction of arbitrary strength. Using a new 
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effective interaction technique we solve in particular the simplest imbalanced case of two spin-up 
and one spin-down fermion (2-|-l) and study its properties. In the so-called fermionized limit of 



^2 infinitely strong repulsive interactions, we solve the problem analytically and find close agreement 

fj with the numerical results. Our findings indicate that the spin-up and spin-down fermions tend to 

phase separate in the trap, which we interpret as a microscopic precursor of ferromagnetism. We 

^ — I 

^ confirm this conclusion by performing calculations for (3-|-l) and (6-|-l) systems for strong repulsive 



interactions. Our predictions are directly addressable in current experiments on ultracold atomic 
few-body systems. 



Introduction. Few-fermion systems are the building blocks of matter. Atoms and nuclei 
are well-known examples, but also systems such as quantum dots, superconducting grains, 
and other nanoscale structures are of great interest . The key to understanding such struc- 
tures is first and foremost the relation between the discrete level structure, due to the finite 
size, and the strength and nature of interparticle interactions. An exciting recent devel- 
opment in atomic physics is the experimental realization of few-body Fermi systems with 
ultracold atoms [U [2] . These setups are extremely versatile as the potential that traps the 
atoms can produce lattices and/or low-dimensional geometries [3], and the atomic interac- 
tion strength may be tuned via the use of Feshbach resonances |1]. The spin-1/2 nature of 
electrons or nucleons is addressable by populating two hyperfine states in the atoms and we 
thus have a direct mapping from the atomic setup to ordinary matter. We will refer to these 
two components as spin up and spin down. 

A seminal contribution of ultracold atomic gas research is the realization of strongly inter- 
acting quantum gases [5HH] using confinement-induced resonances [S] . The Tonks-Girardeau 
(TG) gas [TUHT2] of strongly repulsive bosons that displays fermionic behavior is one such 
example [SHT] . The so-called super-TG (sTG) limit of very strong attractive interactions has 
also been addressed both theoretically [T3H22] and experimentally [8j. Most recently, the 
TG and sTG states have been explored in fermionic systems [231427] . While the two-body 
system in a harmonic trap has a well-known exact solution for any interaction strength, 
known as the Busch model [S^, two-component fermionic few-body systems with more than 
two particles remain largely unexplored. The main difficulty is the handling of very strong 
interactions in the vicinity of the fermionization limit. 

In this letter we face this challenge and consider the experimentally accessible situation of 
a harmonically trapped few-body system in one dimension with N^ spin up and A^^^ spin down 
fermions for the imbalanced case where A^ > A^ = 1. Zero-range interactions of strength 
g are employed between different spin components, while the identical spin particles remain 
non-interacting by the Pauli exclusion principle. This is a few-body analog of the fermionic 
polaron, which is currently under intense study [291 - I3T] . We solve the few-body problem for 
various interaction strengths using a numerical technique inspired by developments in nuclear 
physics [32l |33]. This method shows excellent convergence properties and yields both the 
energy spectrum and full wave functions. In addition, we present an exact solution for the 
A^ = 2 case in the fermionization limit of infinite interaction strength, using an analytical 
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FIG. 1: Energy spectrum for a three-body system (N-^ = 2,N^ = 1) in a harmonic trap as a 
function of the strength, —1/g, of the zero-range interaction between up and down components. 
The cyhnder pfot highhghts the connections of the different states in the non-interacting {g = 0) 
and strong-interacting {\g\ = oo) hmits; on the back and the front of the cyhnder, respectively. 
The states are ordered according to their spatial parity (odd or even). 

model. The classic work of McGuire [3l] solved the untrapped case for arbitrary N^ and g but 
the trapped case has only been solved exactly for N^ = 2 and resonant interaction in three 
dimensions [35j . Our numerical and analytical findings show that around fermionization the 
two spins tend to separate in the trap. We interpret this as a few-body precursor of Stoner 
ferromagnetism [SB] in one dimensional polarized systems. Ferromagnetism is another hotly 
pursued topic at the moment [371 - H2] . As we discuss below, this separation of species should 
be directly observable in current experiments (21 SHI HI] . Other methods have been employed 
recently to similar systems to study energy spectra |l5l US] and also the density profiles |27] . 
Here we present a complementary numerical that convergences extremely fast for multi- 



particle system. We also yield new analytical insights into the problem by obtaining a 
wave function ioi g > that becomes exact at (7 — t- oo and reproduces both densities and 
occupation numbers. 

We concentrate on mostly on the first non-trivial case beyond the two-fermion model, 
which is the N^ = 2, N^ = 1 (2+1) system. Our approach will be described in detail in 
the next paragraph. However, in Fig. [l] we reveal already the numerically obtained energy 
spectrum; plotted on a cylinder where the fermionization limit |(7| — ?■ 00 is on the front 
while weak-coupling \g\ — )■ is on the back. This way of plotting the spectrum emphasizes 
the spectral flow and the connection to the Zeldovich rearrangement effect |17j. The first 
interesting feature is the ground state behavior; it starts as a strongly bound dimer plus a 
particle for g — )■ —00, wraps around the cylinder to the non-interacting limit, \g\ — )■ 0, and 
then becomes energetically degenerate with two other states at (7 —t- 00. One of the states at 
fermionization is a totally antisymmetric combination of all three particles, corresponding to 
placing one particle in each of the three lowest oscillator levels. We will address the physical 
nature of the other states below. 

Methods. We consider a one-dimensional system in a harmonic trap with zero-range 
interparticle interactions for opposite spins. We use harmonic-oscillator (HO) units, in 
which the Hamiltonian is 

^ = E(t + ^)+^ E '^(^^<^-^.^)> (1) 

where lengths are in units of the oscillator length h = yhjmuj., energies in units of the 
trap oscillator energy hio^ cr = ±, and a = —a. In these units, the interaction strength, g, 
becomes dimensionless in units of hub. These units are used for all quantities, unless we 
explicitly state otherwise. The Hamiltonian is parity invariant and one can classify states 
as either even or odd under x — )■ —x. We solve numerically the many-body Schrodinger 
equation in a finite basis constructed from the HO single-particle states |n); each many- 
body basis state is written as \ni..nN^) ® 1^4.) i.e. a product of a HO antisymmetrized state 
of the A'^ spin up particles and a HO single particle state for the spin down particle. The 
corresponding HO energy is (rii + . . . + un + ^4, + y) where N = N^ + Ni. The model space 
can then be defined by an upper limit, i.e. rii + . . . + n^ + n^ < ritot- Since we are only 
interested in the intrinsic dynamics of states, a Lawson projection term is used to push away 
many-body solutions corresponding to excitations of the center of mass motion. 



t. JJ 




1 1 




k 
^ 


•— • gs energy (eff. int.) 


4.30 


^ 
% 
s 


•■ -• gs energy (bare int.) 


4.25 


* 

V 

^ 

^-«. 


- 


"4.20 

4s; 


- -.. 


§4.15 

0) 

i5 4.10 


" " ^ 


4.05 


- 


4.00 




:} Qi; 


1 1 1 1 



20 



30 



40 



50 



60 



n 



tot 



FIG. 2: Evolution of the ground state energy in the 2 + 1 system for 1/g = 0.01 as a function of 
the size ntot of the many-body model-space (excluding the center-of-mass contribution) . The green 
dashed line is the result for the bare Hamiltonian, while the blue solid line is obtained using the 
two-body effective interaction Vp^. In this latter case, ntot = np + 2, where np is the truncation 
of the two-body space. 



Instead of the bare zero-range interaction in ([I]), we consider an effective two-body in- 
teraction in order to speed up the convergence of the eigenenergies with respect to the size 
of the many-body basis. The effective potential Vp^ is constructed in a truncated two-body 
space P, defined as the set of two-body relative HO states whose radial quantum num- 
ber are smaller or equal to a cutoff np. The effective force Vp^ is designed such that its 
solutions correspond to exact solutions given by the Busch formula [28]. Using a unitary 
transformation, we construct a two-body effective Hamiltonian H'^ as [32] : 
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where Epp is the diagonal matrix formed by the np + 1 lowest exact energies given by the 
Busch formula, and Upp is the matrix whose rows are formed by the corresponding eigenvec- 
tors projected on P. The effective interaction Vp^ is obtained from H^ by subtracting the 
HO potential. For each cutoff np, we diagonalize the many-body Schrodinger equation with 
Vp^ and increase ntot until convergence of the many-body energies is reached [32]. We find 
that ntot = rip -|- 2 is sufficient to capture the properties of the effective interaction. With 
this choice, we can then study the energy convergence as a function of ntot- By construc- 
tion, this unitary transformation approach will reproduce exact bare Hamiltonian results, 
both energy spectrum and wave functions, in the np — )■ oo limit. The excellent convergence 
property of our method is demonstrated in Fig. [2] that shows the evolution of the ground 
state energy of the 2-1-1 system as a function of the size of the model space ntot using both 
the bare zero-range potential and the effective interaction Vp^. Results in Fig. obtained 
with Vp^ show a fast convergence of the g.s energy: for ritot = 20, the energy is less than 
0.0005 hu from the converged result. On the contrary, using the bare potential the decrease 
of the g.s. energy as ntot increases is slow and in Fig. |2] convergence is not reached; indeed, 
for the largest model space (ntot = 60) considered here, results are still ~ 0.16 hu above the 
converged value. 

Imbalanced systems. In Fig. |3|we show the 2-1-1 spectrum for odd and even parity states. 
The horizontal lines correspond to totally antisymmetric wave functions, which are non- 
interacting. E.g., the lowest such state, at E = Ahw, corresponds to having one particle in 
each of the three lowest HO states. At (? = oo it becomes degenerate with two interacting 
states. Note also that we have many molecular branches close to (? = oo for (7 < (dimmed 
curves in Fig. pi). Starting from g ^ 0~^ (far left in the figure) and following the odd g.s. we 
see that it makes a jump at g = 00 before becoming non-interacting at (7 — )■ 0^ (far right). 
This is an analog of the repulsive branch for untrapped polarons |30l [31] . However, the jump 
endured by the odd and even states that become degenerate at (7 = 00 is quite different. For 
comparison, we plot the two-body Busch results shifted by the energy of a free spectator 
particle (dashed blue line in Fig. pi), which turns out to be almost identical to the even 
parity state at low energy. This even parity state therefore has an atom-dimer structure, 
with almost no interaction between atom and dimer. This has also been observed in three- 
dimensional traps IH] . A particularly interesting feature of the interacting states, as they 
cross over to the repulsive branch, is their density distributions that we show in Fig. |4] for 
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FIG. 3: Energy spectrum for the 2+1 system (excluding the center-of-mass contribution). Diverg- 
ing states on the attractive side are dimmed. The Busch-model for a 1+1 system plus the energy 
of a spectator particle has been added for comparison. 

the odd ground state. While they are approaching the ferniionization limit \g\ -^ oo for the 
total density, the spin-resolved densities demonstrate a distinct phase separation in the trap. 
This can be seen as a precursor of ferromagnetic behavior in a one-dimensional few-body 
context. This example indicates that the fermionization limit is different for two-component 
fermions as compared to fermionization of bosons. 

Analytical model. We now present an analytical approach that captures the behav- 
ior of the wave function exactly at (7 = 00. We first define Jacobi coordinates, x = 
(xii- — X2-|-)/"\/2 and y = a/2/3x34 — {xi^ + X2^)/VQ, and use these to obtain the spheri- 
cal variables, r = A/aJM-y^ and tancj) = y/x. For g = 0, the (unnormalized) eigenstates 
are r'^L(^(r^)e~'' /^cos(/i0) and r^L'^{r'^)e~^ /^sin(/i0), where /x and u are non-negative in- 
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FIG. 4: Density distributions of the 2+1 ground state as a function of the (repulsive) interaction 
strength, g. Panels (a) and (b) show the spin-separated density distributions for the impurity (spin- 
down) and majority (spin-up) particles, respectively. Panel (c) corresponds to the total density. 
The density of the non-interacting totally antisymmetric state is plotted for comparison together 
with the analytical model for the odd-parity ground state at fermionization. 

tegers and L'^{z) is the associated Laguerre polynomial. The corresponding energies are 
E = 21/ + fi + 1. We need only consider the wave function in the interval < < tt/2 
since one can use the Pauli principle and parity invariance to extend to < < 27r. At 
= 7r/6 opposite spins overlap. The full solution can thus be obtained by matching the 
wave function and its derivative on the line (p = ir/G. We have 



Fi(r, 7r/6) = F2(r, 7r/6) and 
1 9Fi(r,7r/6) 1 dF2{r,Ti/6) 



-Fi(r,7r/6), 



(3) 
(4) 



2r2 d(t) 2r2 dn y/2r 

where Fi and F2 are solutions for < < 7r/6 and tt/6 < < 7r/2 respectively. For 
g ^ these equations are complicated to solve, but by introducing an ad hoc rescaled 
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strength parameter go = gr we decouple the equations and can write Fi{r,(f)) = i?.j(r)\l/j(0) 
for i = 1,2, where Ri{r) = r^L^(r^)e~^' /^. The eigenfunctions and eigenvalue equations 
can now be obtained by using the free angular solutions A cos(/i0) + B sin(/i0) and the 
conditions in Eqs. dsl) and Ml). The nature of the three-fold degeneracy at ferniionization 
seen in Fig. [1] and Fig. [s] comes from the odd and even solutions of these equations |19] , 
while the non- interacting has the structure cos(30). 
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FIG. 5: Occupations numbers as a function of the interaction strength for the ground state of the 
2+1 system. Panels (a) and (b) show the spin-separated occupation numbers, while (c) corresponds 
to the total occupation number. 

Discussion. Current experiments on few-fermion systems [T[ |2] can study the structures 
that we predict by performing tunneling measurements that map out the occupancies of 
the few-body wave function. By varying g one can explore the structure on both sides of 
the resonance [2]. It is possible to go diabatically from the repulsive ground state and onto 
the repulsive branch on the g < side since the overlap with the atom+dimer molecular 
branch is small. It is thus possible to investigate a large part of the parameter space. One 



observable is the occupancies of different single particle levels in the trap, which we show in 
Fig. |5j By selective ejection of the minority particle it is possible to measure the majority 
occupation numbers and a preliminary comparison to experimental data shows agreement 
with our predictions [50] . 

The approach that we have presented can be applied to larger systems. In Fig. |6]we show 
the spin-resolved densities for the ground states of the 3+1 and 6+1 systems as a function 
of the (repulsive) interaction strength. These spin densities show the same phase separation 
behavior in the limit of fermionization as the 2+1 case in Fig. |4j This is in fact found for any 
trapped N+1 system [51]. These general structures can also be investigated by tunneling 
0152]. 
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FIG. 6: Spin-resolved densities for the 3+1 and 6+1 systems, cf. Fig. |4| Panels (a) and (c) show 
the distribution of the impurity particle, while panels (b) and (d) show the majority. 
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